High order elastic terms and general paradigm of the elastic interaction between 
colloidal particles in the nematic liquid crystals. 
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Theoretical description of the elastic interaction between colloidal particles in NLC with incor- 
poration of the higher order elastic terms beyond the limit of dipole and qudrupole interactions is 
proposed. The expression for the elastic interaction potential between axially symmetric colloidal 
particles, taking into account of the high order elastic terms, is obtained. The general paradigm 
of the elastic interaction between colloidal particles in NLC is proposed so that every particle with 
strong anchoring and radius a has three zones surrounding itself. The first zone for a < r J$ 1.3a is 
the zone of topological defects; the second zone at the approximate distance range 1.3a J$ r 5; 4a is 
the zone where crossover from topological defects to the main multipole moment takes place. The 
higher order elastic terms are essential nere (from 10% to 60% of the total deformation). The third 
zone is the zone of the main multipole moment, where higher order terms make a contribution of 
less than 10%. This zone extends to distances r 4a = 2D. 

The case of spherical particles with planar anchoring conditions and boojums at the poles is 
considered as an example. It is found that boojums can be described analitically via multipole 
expansion with accuracy up to 1/r 7 and the whole spherical particle can be effectively considered as 
the multipole of the order 6 with multipolarity equal 2 6 = 64. The correspondent elastic interaction 
with higher order elastic terms gives the angle 9 m i n = 34.5° of minimum energy between two contact 
beads which is close to the experimental value of 8 m i n = 30°. In addition, high order elastic terms 
make the effective power of the repulsive potential to be non-integer at the range 4.5 < 7 e // < 5 for 
different distances. The incorporation of the high order elastic terms in the confined NLC produce 
results that agree with experimental data as well. 



PACS numbers: 



I. INTRODUCTION 



Anisotropic properties of the nematic liquid crys- 
tals (NLC) give rise to a new class of colloidal elas- 
tic anisotropic interactions that never occur in isotropic 
hosts and result in different structures of colloidal par- 
ticles: linear chains [l], H|, inclined chains with respect 
to the director and quasi 2D nematic colloids [||- 

[l3j. Theoretical understanding of the matter in the 
bulk NLC is based on the multipole expansion of the 
director field and has deep electrostatic analogies. Untill 
now, all theoretical models dealt with only the first three 
terms in multipole expansion: Coulomb- like [20j . dipole 
and quadrupole fl7j-[26l]. Almost all experiments are 
made with axially symmetric colloidal particles (primar- 
ily spherical) which carry only dipole and/or quadrupole 
clastic moments. But considering only these two terms 
cannot explain quantitatively any of the observed struc- 
tures. For instance the droplets with tangential boundary 
conditions make an angle of 30° with the alignment axis 
of the liquid crystal , which is along the vertical 
axis. However, the quadrupole interaction gives the an- 
gle, for which long-range attraction is maximized to be 
approximately 49°. Therefore the origin of the existing 
structures must be ascribed to the short-range effects, 
not explicitly included in the theory. 
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In the current paper it is found that high order multi- 
pole terms play a very important role in the short-range 
effects and in the formation of colloidal structures. Ac- 
tually we find that there are three zones around each 
colloidal particle: the first zone is the zone of topological 
defects where non-linear terms are essential. It has the 
approximate size of 0.2a — 0.3a of the particle radius a, 
so that it is concentrated on the distances a J$ r J$ 1.3a. 
The second zone is the intermediate zone, where all pos- 
sible from the symmetry point of view clastic terms, 
are born simultaneously and higher order elastic terms 
are essential (from 10% to 60% of the total deforma- 
tions). It is concentrated at the approximate distances 
1.3a ;$ r J$ 4a. And the third zone is the zone of the 
last multipole moment, where higher order terms make 
contribution less than 10% and only the last multipole 
moment has the dominant value. This third zone extends 
to distances r Si 4a = 2D. 

We will now consider the case of spherical particles 
with planar anchoring conditions as an example. Such 
particles have topological defects called boojums at the 
poles. We find that boojums can be effectively described 
via multipole expansion with accuracy up to 1/r 7 and 
the whole spherical particle can be effectively consid- 
ered as multipole of the order 6 with multipolarity equal 
2 6 = 64. The correspondent elastic interaction between 
two beads with higher order elastic terms gives the an- 
gle B m in = 34.5° of minimum energy between two con- 
tact beads which is close to the experimental value of 
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II. INCORPORATION OF THE HIGHER 
ORDER ELASTIC TERMS INTO THE THEORY 

Let's now consider axially symmetric particle of the 
micron or sub-micron size which may carry topological 
defects such as hyperbolic hedgehog, disclination ring 
or boojums. In the absence of the particle the non- 
deformed state of NLC is the orientation of the director 
n||z, n = (0,0, 1). The immersed particle induces defor- 
mations of the director in the perpendicular directions 
rip, /i = x,y and make director field n w (n x , n y , 1). The 
bulk energy of deformation may be approximately writ- 
ten in the harmonic form: 




with Euler-Lagrangc equations of Laplace type: 

= (2) 

Then the director field outside the particle in the infinite 
LC has the form n x (r) = p^ + 3cff , n y (r) = p± + 3c*f 
in the simplest case with p and c being dipole and 
quadrupole elastic moments. The anharmonic correc- 
tion to the bulk energy is F an h ar = ^ J d 3 x(Vn z ) 2 sa 
Y J d 3 x(X7n 2 L ) 2 which changes EL equations to be: 

An M + -n^Anl = (3) 

If the leading contribution to is the dipolar term then 
anharmoic corrections are of the form r^/r 7 and high or- 
der terms of the order up to 1/ r 5 can effectively influence 
on the short-range behaviour and should be equally con- 
sidered. The same if the leading contribution to n M is the 
quadrupolar term then anharmonic corrections are of the 
form r^/r 10 and high order terms of the order up to l/r 8 
can effectively influence the short-range behaviour. 

In the general case, the solution of the Laplace equa- 
tion for axially symmetric particles has the form: 

N 1 
n M = X;oi(-l)'^ _1 - (4) 

where a/ is the multipole moment of the order I and 
2 l is the multipolarity; A - is the maximum possible 
order without anharmonic corrections. For the dipole 
particle N = 4, for the quadrupole particle N = 7. 
So ai = p is the dipole moment, a 2 = c - is the 
quadrupole moment. Actually, all odd coefficients are 
equal to zero for quadrupole particles 03 = 05 = a? = 
because of the horisontal symmetry plane so that it can 
be limited with N = 6. All nonzero coefficiants a/ 
are unknown quantities. They can be found as asymp- 
totics from exact solutions or from variational ansatzes. 
Strictly speaking these coefficients are functions of the 
anchoring coefficient W and surface elastic constants 
ai = ai(Wa/K, K24/K, Ki 3 /K) as well as the Frank elas- 
tic constant K and particle radius a. But we don't ex- 
actly know this dependence. From the other side the 




FIG. 1: (Color online) Horizontal projection of the director 
n p = n± at the surface of the spherical particle r = a with 
planar anchoring conditions for < 9 < 90° (see Fig[2]d). 
According to the solution (|10[) . Blue line 1 corresponds to 
(b 2 ,&4,&e) = (-0.36,-0.023,-0.00018), blue dashed line 2 
corresponds to ba — 0, (62, &4, b§) — (—0.32, —0.007, 0) and red 
line 3 corresponds to the pure quadrupole term for b±, ba = 0, 
(b 2 ,b 4 ,b 6 ) = (-0.28,0,0). 



coefficients a/ may be found as fitting parameters in the 
interaction potentials for each particular case. 

In order to find the energy of the system: particle(s) 
+ LC , it is necessary to introduce some effective free 
energy functional F e ff so that it's Euler-Lagrange equa- 
tions would have the above solutions (g]) . In the one con- 
stant approximation with Frank constant K the effective 
functional has the form: 

Feff J ^jiZ^^^^AKx)^ 1 ^! 

(5) 

which brings Euler-Lagrange equations: 

N 

Atv = 4 7 r^(-l)'- 1 9 Al 9i- 1 A ; (x) (6) 
;=i 

where A;(x) are multipole moment densities, fi = x,y 
and repeated y, means summation on x and y like d^n^ = 
d x n x + d y n y . For the infinite space the solution has the 
known form: 




^(-lyc^-^x') (?) 



If we consider Ai (x) = a/<5(x) this really brings solution 
(j3J). This means that effective functional ([S]) correctly 
describes the interaction between the particle and LC. 

Consider N p particles in the NLC, so that ^4;(x) = 
J^a^x — Xj), i = 1 -f- N p . Then substitution ([7]) 
into F eff © brings: F eff = U self + u interaction where 
jjseif _ . jj^ el f t here U- el ^ is the divergent self energy. 
Interaction energy u interaction = J2i<j U™ 1 ■ Here Uff 
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FIG. 2: (a) Taken from 0. Colloidal particle (2a = 4.5/zm) 
immersed in nematic liquid crystal in a planar cell as seen 
under a polarizing microscope. Two boojums at the poles of 
the sphere confirm tangential alignment of the director at the 
sphere surface, (b) Taken from Aggregation of particles 
was observed at an angle of approximately 6 ~ 25° — 35° with 
respect to the average director orientation no||2.(c) Intensity 
profile of the nematic field obtained by numerical calcula- 
tion shown by grayscale plots of n 2 p from the solution (| 10|) for 
(b 2 ,b 4 ,b 6 ) = (-0.36, -0.023, -0.00018). In the black regions 
the director aligns along the z axis, while deformations n 2 p are 
maximal in the white regions, (d) Chart of the director field 
for the boojums configuration. 



is the elastic interaction energy between i and j particles 
in the unlimited NLC: 

Uif = 4nK ]T a,a{,(-l)' ^±P l+l ,(cos9) (8) 
i,i'=i 

Here unprimcd quantities ai are used for particle i and 
primed a[, for particle j, r = |x; — Xj|, 9 is the angle 
between r and z and we used the relation Pi(cos9) = 
(—l) lL ^-d l z ^ for Legendre polynomials Pi. It is the gen- 
eral expression for the elastic interaction potential be- 
tween axially symmetric colloidal particles in the un- 
bounding NLC with taking into account of the high order 
clastic terms. 

The case of confined NLC means just replacement of 
r = | x--x- | w ith the Green's function G(x, x') (see 0, 
l24j ]) , which satisfies equation A x G(x, x') = — 47r<5(x— x') 
for x, x' G V (V is the volume of the bulk NLC) and 
G(x, s) = for any s of the bounding surfaces S. Then 



formula §E§ for the confined NLC has the form: 

N 

jjint, con fined _ _^ JT j i ' £ a^d^d^d 1 / "^(x,, x'j) 
l,l' = l 

(9) 

For dipolc particles the sum is limited to nonzero 
ffli) G2, a-3 and a^. For quadrupole particles (beads with 
boojums and Saturn ring configuration) the sum is lim- 
ited to nonzero 02,04 and a^. All coefficients may be 
presented as a/ = 6;a i+1 with a being the radius of the 
particle and bi are just dimensionlcss parameters. 

Below we consider a spherical particle with planar an- 
choring conditions at the surface as an example. Then 
the director field ((4]) can be presented as n a = a2d n d z ^ + 
a^dad^ + a 6 d p d z ^. Let's introduce dimensionless dis- 
tance r =>■ ar, then the horizontal projection n p = n±_ 
(?ip = n u n u ) of the director has the form: 

, 3sin9cos9 , 105sin9cos 3 9 — A5sin9cos9 



, I0395sin9cos 5 9 - 9450sin9cos 3 9 + 1575sin9cos9 
+ 6 6 ~ 7 

(10) 

where 9 is the angle between z and r, n x = n p cos(p, 
ii y = Upsintp and ip is the angle between p and x. It 
is obvious that &2 < for the boojums configuration. 
In the paper [24| it was found that 62 ~ —0.28 for the 
experiment and we have two unknown variation pa- 
rameters 64 and &6- Physical limitations for these coeffi- 
cients may be formulated in the following way: \n p \ < 1 
for all distances r > 1; n p (r,9) should have only one 
minimum/maximum as a function of 9 for < 9 < ir/2 
and for all r > 1 and the correspondent energy of in- 
teraction ([5]) should agree with all known experimental 
results as much as possible. We found that all these 
conditions are satisfied in the best way for the values 
O2, 64, b 6 ) = (-0.36, -0.023, -0.00018)(scc below). 

The dependence n p {9) on the spherical surface r = 1 
is depicted on the FigQ] for the values (62,^4,^6) = 
(-0.36, -0.021, -0.00011) (blue line 1). Topological de- 
fects called boojums are located at the poles of the parti- 
cle (see Fig(2]a). The intensity profile of the nematic field 
obtained by numerical calculation shown by grayscale 
plots of n 2 p from the solution (|10p is presented on the 
FigEJc . In the black regions the director aligns along 
the z axis, while deformations n 2 p are maximal in the 
white regions. We see that pictures FigJ5]a and Fig(21c 
are quite similar. This means that solution (|10j) gives the 
correct analitical description of the boojums near the sur- 
face of the spherical particles up to the 1/r 7 order. Of 
course there will be some corrections to the solution from 
the anharmonic term, but they will decrease faster than 
1/r 7 as it is seen from the equation ^ . 

We will discuss, in more detail, the limits of applica- 
bility of the higher order elastic terms as well as make a 
more profound estimation of nonlinear terms in the Sec 
III. 
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FIG. 3: (Color online) Experimental values (taken from 0]) 
of interparticle repulsion force F between two beads of di- 
ameter D = 2a = 4.4/im as a function of rescaled dis- 
tance x/D. Blue line 1 is calculated as F = — U' r from 
JTIJ for 6 = tt/2, K = 7pN (5CB) and (6 2 ,&4,6 6 ) = 
(-0.36, -0.023, -0.00018), blue dashed line 2 corresponds to 
be = 0, (& 2 ,&4,6 6 ) = (-0.32,-0.007,0) and red line 3 cor- 
responds to the pure quadrupole-quadrupole interaction for 
& 4 , b 6 = 0, (6 2 , b 4 , b 6 ) = (-0.28, 0, 0) . 



A. The effect of the high order terms on the 
angular dependence of the interaction potential 

The correspondent energy (J8]) of elastic interaction be- 
tween two spheres with boojums has the form (all coeffi- 
cients bi = besides (62, 64, b$) and N = 6): 



U 



AnaK 



-(2& 2 & 6 



Pa(cosL 
P%{cos9] 



26-2^46! 



26 4 6 6 10! 



Pe {cosl 
P w {cot 



(11; 
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FIG. 4: (Color online) Blue line 1 corresponds to the angular 
dependence of the interaction potential in kT units for 
two close beads at the distance r — 2a, a = 2.2/im, K = 7pN 
for (62,64,60) = (-0.36, -0.023, -0.00018). The potential has 
minimum at 9 = 34.5°. Blue dashed line 2 corresponds to the 
case 6 6 = 0, (6 2 ,6 4 ,6 6 ) = (-0.32,-0.007,0). The potential 
has minimum at 6 — 44°. Red line 3 corresponds to the pure 
quadrupole-quadrupole interaction for 64, be = 0, (62, 64, be) = 
(-0.28, 0, 0) with minimum at 9 = 49° . 
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FIG. 5: (Color online) The angular dependence of the in- 
teraction potential in kT units for two beads at differ- 
ent distances. Broun dashed line 1 corresponds to r = 2.5a, 
black dotted line 2 for r = 3a, black line 3 for r = 4a, red 
line 4 for r = 5a. Here radius a = 2.2fim, K = 7pN and 
(6 2 , 64, b 6 ) = (-0.36, -0.023, -0.00018). 




where r is measured in radius a units and r > 2. The 
correspondent force F = — U' r of repulsion between two 
beads of diameter D = 2a = iAfxm in 5CB (K = 7pN) 
for 9 = tt/2 is plotted on the FigJ3] We see that 
experimental values of the repulsion force may be fit- 
ted with three different sets of parameters: blue line 1 
corresponds to (6 2 ,6 4 ,6 6 ) = (-0.36,-0.023,-0.00018); 
blue dashed line 2 corresponds to the case b§ = 0, 
(62AA) = (-0.32,-0.007,0). Red line 3 corre- 
sponds to the pure quadrupole-quadrupole interaction 
for 64, b 6 = 0, (6 2 ,6 4 ,6 6 ) = (-0.28,0,0). It is very in- 
teresting that all these three set of parameters fit the 
data very well on the distances 1 < r/D < 1.6. Actu- 
ally the experimental values on the Fig|3] were found in 



the homeotropic cell with width L = 1.8D = 7.9fj,m Q. 
It was found in 0, [22|, [23| that confining effects become 
essential for distances of more than r > 0.9L so that we 
can use approximation of the unbounding NLC (|11[) for 
the distances r < 0.9L = 1.6D (for L = 1.8D). 

Despite the fact that three different set of parame- 
ters give almost the same values of the repulsion force 
in the perpendicular direction 9 = tt/2, they produce 
very different pictures for the angular dependence of the 
interaction potential. The angular dependences of the 
interaction potential for two close contact beads at the 
distance r — 2 is depicted on the FigUJ It is clearly 
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FIG. 6: (Color online) Dependence of the minimum en- 
ergy angle O m in on the rescaled distance r/a between 
two beads for the potential Hll[) with (62,64,6s) = 
(-0.36, -0.023, -0.00018) (blue line 1), minimum energy an- 
gle 9 = 49° for pure quadrupole interaction (red line 2). 



Jeff 



2 




FIG. 7: (Color online) Dependence of the effective power 
on the rescaled distance r/a between two beads in the per- 
pendicular direction 8 = tt/2 for the potential with 
(62, 64, &e) = (-0.36, -0.023, -0.00018). 



seen that two beads with set of parameters (62, 64, 6g) = 
(-0.36,-0.023,-0.00018) (FiggJ blue line 1) produce 
the interaction potential (jllj) which has the minimum 
at the angle = 34.5° that is very close to the results 
observed earlier in experiments [3, 4] (sec FigJ^Jb). The 
set of parameters (6 2 , 64, &e) = (-0.32, -0.007, 0) ( FigH 
blue dashed line 2) produces the potential with the mini- 
mum at the angle 9 = 44° . Red line 3 corresponds to the 
pure quadrupolc-quadrupole interaction for 64,66 = 0, 
(62,64,6s) = (—0.28,0,0) and the correspondent interac- 
tion potential (jllj) has minimum at 9 = 49°. 

The angular dependence of the elastic interaction po- 
tential for different distances r is depicted on the FigJS] 
The minimum energy angle increases from 9 = 34.5° to 
9 = 48° with increase of the distance from r = 2a to 
r = 8a that is shown on the Fig|5] This corresponds 
to the experimental results of Q where minimum energy 
angle 9 was found to be changed from 9 = 30° to 9 = 48° 
with increase of the distance between particles. 

So we come to the conclusion that high order elastic 
terms have very profound influence on the angular depen- 
dence of the interaction potentials at the short distances 
between particles which agree with experimental results. 



B. The effect of the high order terms on the 
effective power 

Simple electrostatic analogy developed in [lTj predicts 
that elastic forces are proportional to F oc r~ 4 ,r -5 ,r -6 
for different types of clastic interactions. But many ex- 
periments give non-enteger power dependence F oc r~ s 
with S = 3.6 in [HI, S = 4.6 in for different di- 

rector configurations. In the paper [161 ] this descrepancy 
was succesfully fitted with help of possible contribution 
of higher order terms in multipolc expansion of F. We 



argue, as well, that high order elastic terms make an ef- 
fective power to be non-integer in the range of severel 
percents. 

Let's consider the potential (fTTj) for 9 = 
7r/2 and the set of parameters (62,64,65) = 
(-0.36,-0.023,-0.00018). This potential is repul- 
sive elsewhere. Let us present it, approximately, in the 
form of power law dependence with some effective power 
that depends on the distance, i.e.: 

C 

rrwall w (1 p\ 

^ same Mom ~ V / 

where -y eff may be found as 7 e// = -f^§7 = -U' r jj. 
Fig 111] shows dependence of such effective power on the 
dimensionless distance r/a. We see that on small dis- 
tances 2 < r/a < 3 effective power j e ff decreases from 
5.2 to 4.5 and then it increases from 4.5 to 5 for r/a > 3. 

C. The effect of the high order terms in the 
confined cell 

The formula (O gives the elastic interaction poten- 
tial in the confined NLC. Lets's consider homeotropic 
ncmatic cell with thickness L. The Green function 
G^^(x, x') coincides with the Green function G(x, xQ of 
the two conducting walls in the electrostatics (see |27j): 



n— 1 m— — oo 

. mrz' ,mr P< ,nir P< 
xsm— — J m ( - )K m { - ) (13) 

Here heights z, z', horizontal projections p<,p> and 
I m ,K m are modified Bessel functions. Then using of 
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der to do this we need to estimate anharmonic energy 
term F anhar = f / <Px(Vn z f = f J d 3 x&^j = 
J d 3 x fanhar (x) and compare it with the harmonic term 
Fhar = y J d 3 x(Vn u ) 2 = J d 3 xf har (x-)- In addition, we 
need to compare the contribution of the high order elas- 
tic terms to the total deformation and compare it with 
the contribution of the main multipole term. 

If we substitute the solution (TIU]) with (62,64,65) = 
(-0.36, -0.023, -0.00018) into F anhar and F har we re- 
ceive after numerical integration F nar = 3.239i^a, 
Fanhar = 1.065-Ka so that total deformation energy is 
Fdeform = 4.3-Ka. Let's analize where the anharmonic 
energy term is the most localized. To do this we intro- 
duce the ratio of the anharmonic and harmonic energy 
densities: 



x/D 



FIG. 8: (Color online) (Color online) Experimental data taken 
from 0] - energy of elastic interaction between two spherical 
particles with diameter D — 2a = 4.4/im in the homeotropic 
cell with thicknesses L — h — 6.5/im and h — 8fim. The solid 
blue line 1 and green line 2 are calculated according to the 
formula (|14[) for z = z' = L/2. The parameters (62,64,66) = 
(-0.36, -0.023, -0.00018) fit both thicknesses pretty well in 
the energy scale 1-1000A:T. 



(|9"| brings the elastic interaction between two beads with 
boojums in the homeotropic cell : 



t Thorn . , 

16naK \LJ ^ ' L L L 
x \b\ + 2b 2 b i (mr) 2 (£) 2 + (bl + 2b 2 b 6 )(mr) A 

+ 2bMnn) 6 (|) 6 + blinvf [ff] (14) 

with p being the horizontal projection of the distance 
between the particles. 

Fig J5] demonstrates the application of this formula (fT4")) 
for the repulsion potential between two spherical par- 
ticles (with planar anchoring on the surface providing 
quadrupole director configuration) with diameter D = 
2a = 4Afim in the center of homeotropic cell (z = z' = 
L/2) with thicknesses L = 6.5utti and L = 8/xm ( exper- 
imental data are taken from [7| )■ It is seen that the set 
of parameters (62,64, 6 6 ) = (-0.36, -0.023, -0.00018) fit 
both thiknesses pretty well in the energy scale l-1000fcT. 



£ anharrn 



(x) 



fanhar (^) 

jw(x) 



4(Vrv)2(l-ni) 



(15) 



The Fig IH1 demonstrates e an harm{z, %) in the plane ZX 
(no||z). It is clearly seen that s an harm(z, x) is localized 
in the area of 0.2a — 0.3a from the particle surface and 
it becomes rapidly less than 10% further (see FigfTU] as 
well. ) 

Fia llll demonstrates perpendicular projection of the 
director field n p = n± at the different distances from 
the center of the particle according to the solution (flU)) 
with (62,64,65) = (-0.36,-0.023,-0.00018). It is large 
enough 0.4 < n p < 0.8 on the distances 1 < r/a < 1.2 
and non-linear corrections to the solution f| 10[) may be es- 



sential here. But further it becomes smaller n p < 0.3 for 
1.3 < r/a so that condition of harmonic approximation 
n u 1 is satisfied and the director is very well described 
by the solution (|10[) . Thus high order terms are born 
and are essential for distances 1.3 < r/a on the left side. 
Let's estimate where the end of their influence is on the 
right side. 

Let up 10 - be the part of the director deformation pro- 
duces by the high order elastic terms (see (fTU)) ): 



h.o. 



105 sin8 cos 9 — AbsinOcosO 



10395sm6»cos 5 6» - 9450sm6»cos 3 



1575sin9cos9 



+6, 

r 

(16) 

Let's introduce the fraction rjh.o. of the deformation 
produced by the high order elastic terms in the total 
deformation: 



III. THE INFLUENCE OF NONLINEAR 
TERMS. THE LIMITS OF APPLICABILITY OF 
THE HIGHER ORDER ELASTIC TERMS. 

In this section we want to discuss the limits of ap- 
plicability of the higher order elastic terms. In or- 



Vh-aX*) = 



7^°-(x) 



n p (x) 



(17) 



where the total deformation n p (x) is defined in (TIT))) . 

Figfl2l demonstrates the fraction r)h. X r ) f° r two dif- 
ferent directions: 8 = (blue line 1) and 9 = ir/2 (brown 
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FIG. 9: (Color online) The ratio of the anharmonic and har- 
monic energy densities E an harm(x) (|15p in ZX plane around 
the spherical particle of the unit radius. 



anharm 




20 40 60 80 qO 

FIG. 10: (Color online)The ratio of the anharmonic and har- 
monic energy densities £ a nharm{9, r) in the first quarter of the 
FigO for different distances from the center of the spherical 
particle (8 is the angle between r and z). Blue line 1 corre- 
spond to the distance - r/a = 1.06, lilac line 2 - r/a — 1.1, 
brown line 3 - r/a — 1.2, green line 4 - r/a = 1.3 



line 2) with (b 2 , h, b 6 ) = (-0.36, -0.023, -0.00018). It is 
clearly seen that this fraction is large enough (from 10% 
to 60%) in the range r/a < 4 and it is less than 10% for 
r/a > 4 where the main multipole term (quadrupole in 
this case) plays the dominant role. Therefore we can say 
that high order elastic terms play an important role in 
the range 1.3 < r/a < 4. 



IV. GENERAL PARADIGM OF THE ELASTIC 
INTERACTIONS BETWEEN COLLOIDAL 
PARTICLES IN NLC 

The results obtained above help us to formulate the 
following picture or paradigm of the elastic interaction 




FIG. 11: (Color online) Horizontal projection of the direc- 
tor rip = n± at the different distances from the center of the 
particle. According to the solution (|10|) with Q)2,b4,be) = 
(-0.36, -0.023, -0.00018). Blue line 1 correspond to the dis- 
tance - r/a = 1, lilac line 2 - r/a — 1.1, brown line 3 - 
r/a = 1.2, green line 4 - r/a — 1.3 




FIG. 12: (Color online) The fraction of the high order elastic 
deformations rjh.o.( r ) C3) m the total deformations around 
the spherical colloidal particle. Blue line 1 corresponds to the 
direction 9 — 0, brown line 2 corresponds to the direction 
9 = tt/2 



between colloidal particles in NLC (see FigQ~3]). 

There are three different zones around each colloidal 
particle. The first zone is the zone of topological defects 
(brown zone 1 on the Fig |13[) . Non-linear terms are very 
essential in this zone, the EL equation is non-linear. So 
that the principle of superposition does not work in the 
first zone. The size of the first zone is about 0.2a — 0.3a 
of the particle radius a, so that it is concentrated at the 
distances a < r ;$ 1.3a. This is in line with results ob- 
tained for the hyperbolic hedgeho g an d Saturn ring di- 
rector configurations. In the paper [17( it was found that 
hyperbolic hedgehog is located at the distance r = 1.22a 
and Saturn ring is located at the distance r = 1.08a. 

The second zone appears just after the first zone (dark 



8 




FIG. 13: (Color online) General structure of the nematic elas- 
tic field around the colloidal particle. Brown zone 1 - is the 
zone of topological defects and anharmonic deformations. The 
principle of superposition does not work here. Dark green 
zone 2 - is the zone of harmonic elastic terms of all possible 
orders. High order elastic terms are essential here (from 10% 
to 60% of the total deformation). Light green zone 3 - is the 
zone of the main elastic term. 



green zone 2 on the Fig |13[) . In this zone anharmonic 
terms vanish and harmonic elastic terms of all possible 
from the symmetry point of view orders are born simul- 
taneously. The director field here can be presented as 
multipolc expansion of all possible orders and the prin- 
ciple of superposition is valid. All clastic terms coexist 
here and higher order elastic terms are essential (from 
10% to 60% of the total deformation). It is concentrated 
at the approximate distance range 1.3a ^ r ;$ 4a. The 
second zone is the zone where crossover from topological 
defects to the main multipole moment takes place. 

And the third zone is the zone of the main multipole 
moment, where higher order terms make contribution less 
than 10% and only the first nonzero multipole moment 
has the dominant value. This zone extends for distances 
r <; 4a = 2D (light green zone 3 on the FigH3|). The 
influence of the high order terms still exists in the third 
zone and has a contribution of approximately several per- 
cent. For instance, it makes the effective power to be 
non-integer like it was shown above. 

V. CONCLUSION 

To conclude, theoretical description of the elastic inter- 
action between colloidal particles with the incorporation 



of the higher order elastic terms is proposed. 

The general paradigm of the elastic interaction be- 
tween colloidal particles in NLC is proposed. Each parti- 
cle has three zones around itself; the first zone is the zone 
of topological defects where anharmonic terms are essen- 
tial and the principle of superposition does not work. 
This zone has the size of about 0.2a — 0.3a of the parti- 
cle radius a, so that it is concentrated at the distances 
a < r ^ 1.3a. 

The second zone appears just after the first zone. In 
this zone anharmonic terms quickly vanish and harmonic 
elastic terms of all possible from the symmetry point of 
view orders are born simultaneously. The director field 
here can be presented as multipole expansion of all possi- 
ble orders and the principle of superposition is valid. The 
higher order elastic terms are essential here (from 10% 
to 60% of the total deformation) and this zone is concen- 
trated at the approximate distance range 1.3a ~ r ;$ 4a. 
It is the zone where crossover from topological defects to 
the main multipole moment takes place. 

The last third zone is the zone of the main multipole 
moment, where higher order terms make a contribution 
of less than 10% and only the first nonzero multipole 
moment has the dominant value. This zone extends for 
distances r 4a = 2D. 

Of course all three zones exist only for particles with 
strong anchoring conditions at the particle surface. The 
first zone is absent for particles with weak anchoring and 
the second zone starts just from the particle's surface in 
this case. 

We consider the case of spherical particles with planar 
anchoring conditions as an example. Such particles have 
topological defects called boojums at the poles. We find 
that boojums can be described analitically via multipole 
expansion with accuracy up to 1 /r 7 and the whole spher- 
ical particle can be effectively considered as the multipole 
of the order 6 with multipolarity equal 2 6 = 64. The cor- 
respondent elastic interaction between two beads with 
higher order elastic terms gives the angle Q m in = 34.5° 
of minimum energy between two contact beads which 
is close to the experimental value of 8 m i n = 30°. As 
well high order elastic terms make the effective power 
of the repulsive potential to be non-integer at the range 
4.5 < 7e// < 5 for different distances. The incorporation 
of the high order elastic terms in the confined NLC pro- 
duce results that agree with experimental data as well. 

The application of the higher order terms for hyper- 
bolic hedgehog and Saturn ring director configurations is 
under the way. The author is grateful to Prof. B.I. Lev 
for fruitful discussions. 
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